Intra-molecular coupling as a mechanism for a liquid-liquid phase transition 
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We study a model for water with a tunable intra-molecular interaction J a , using mean field 
theory and off-lattice Monte Carlo simulations. For all J CT > 0, the model displays a temperature 
of maximum density. For a finite intra-molecular interaction J a > 0, our calculations support the 
presence of a liquid-liquid phase transition with a possible liquid-liquid critical point for water, likely 
pre-empted by inevitable freezing. For J — the liquid-liquid critical point disappears at T — 0. 
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I. INTRODUCTION 

Water has an anomalous density decrease for isobaric 
cooling below a temperature of maximum density (TMD) 
0. Other thermodynamic anomalies — such as the rapid 
increase of the response functions — can be fit by power 
laws with apparent singularity well below the freezing 
temperature Q. Several interpretations for this behav- 
ior have been proposed, but it is unclear which describes 
water or if any describes other anomalous liquids includ- 
ing, among the others, S, Se, Te, Cs, Si, Ge, L C. P, Si02, 

BeF 2 gig! | BUS n 0, n B H 0. 

One of the interpretations, the stability-limit conjec- 
ture [^9), assumes that the limits of stability of the su- 
perheated, supercooled and stretched liquid form a single 
retracing spinodal line in the pressure-temperature (P- 
T) plane. This scenario predicts a divergence of the re- 
sponse functions at the supercooled liquid-to-liquid spin- 
odal M. 



The singularity-free interpretation 21, 22 envisages 



the experimental data represent apparent singularities, 
due to anticorrelated fluctuations of volume and entropy. 
In this scenario, these fluctuations are responsible for the 
TMD line. 

The liquid-liquid phase transition hypothesis |p3f pro- 
poses the presence of a first order line of phase transitions 
separating two liquid phases differing in density, the high 
density liquid (HDL) and the low density liquid (LDL). 
In this scenario the HDL-LDL phase transition, possibly 
ending in a liquid-liquid critical point, is responsible for 
the anomalies. 

Although Refs. || ||], by tuning parameters of the cor- 
responding models, predict smooth transitions from the 
different scenarios, to help elucidate which is the most 
reasonable description for water, we consider a model 
fluid with inter-molecular and intra-molecular interac- 
tions. This model, in the particular case of a zero intra- 



molecular interaction, recovers the model introduced by 
Sastry et al. [Q that predicts the singularity-free sce- 
nario. Our aim is to understand how the presence of an 
intra-molecular interaction changes this prediction. 

We perform analytic calculations in a mean field ap- 
proximation and an off-lattice Monte Carlo (MC) simula- 
tion. Our results show that a non-zero intra-molecular in- 
teraction gives rise to a HDL-LDL phase transition, with 
a possible critical point p3[ , with the liquid-liquid criti- 
cal temperature decreasing to zero and vanishing with the 
intra-molecular interaction. Therefore, at least for this 
model, the singularity-free scenario is obtained only in 
the particular case of a zero intra-molecular interaction, 
while for a finite intra-molecular interaction the HDL- 
LDL phase transition is predicted. General considera- 
tions suggest that the liquid-liquid phase transition for 
water could occur below the glass temperature, i.e. out- 
side the accessible experimental range. 

The paper is organized as follows: in Sec. II we define 
the model defined on a lattice. In Sec. Ill we describe 
the equation of state approach in the mean field approxi- 
mation and we present the mean field results. In Sec. IV 
we introduce the off-lattice model, we describe the MC 
approach and show the simulation results. In Sec. V we 
discuss our results and we give the conclusions. 



II. THE LATTICE MODEL 

The fluid is represented by partitioning the system into 
N cells of equal size. A variable rii is associated with each 
cell i — 1, . . . , N, with n, = 1 if the cell is occupied by a 
molecule, rij = otherwise. 

The inter- molecular interaction |23], 



n 



-e riiUj 

(id) 



J 



Uin i l 

(id) 



(1) 



'Present address: SMC-INFM, Dipartimento di Fisica, Universita 
"La Sapienza", P.le A. Moro 2, 1-00185 Roma, Italy 



has a first term describing the van der Waals attraction 
between molecules, where e > is the energy gain for 
two nearest neighbor (nn) occupied cells and the sum is 
over all the possible nn cells. 

The second term in Eq. (pi) accounts for the dynamic 
network of hydrogen bonds (HBs) formed by liquid wa- 
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ter, with each molecule typically bonded to four other 
molecules at low T |J, with an energy-gain J > per 
HB. We consider cells with size of a water molecule and 
with four arms, one per possible HB. For the molecule 
in the cell i, the orientation of the arm facing the cell j 
is represented by a Potts variable a%j — 1, . . . , q, with a 
finite number q of possible orientations. Two molecules 
in nn cells form a HB only if they are correctly oriented 



i.e., by assumption |22fl, if S a 



1 (*«,. 



1 if 



a = b and S a ^ = otherwise). 

The experimental oxygen-oxygen correlation function 
shows that a HB is formed if and only if the inter- 
molecular distance is within a characteristic range jl). 
Hence, we assume that the formation of a HB leads to a 
local volume expansion |22|, 



V = V + N HB v H B , 
where Vq is the volume of the liquid with no HBs, 



Nhb 



riinjSo 



(2) 



(3) 



is the total number of HBs in the system, and vhb is the 
specific volume per HB. 

Experiments show that the relative orientations of the 
arms of a water molecule are correlated, with the average 
H-O-H angle equal to 104.45° in an isolated molecule, 
104.474° in the gas and 106° in the high-T liquid ||, 
suggesting an intra-molecular interaction between the 
arms. This interaction must be finite, because the an- 

t changes with T, consistent with ab-initio calculations 
and molecular dynamics simulation s [p6| . Hence, we 
introduce the intra-molecular term 
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where for each of the 4 C2 = 6 different pairs (k, l)i of the 
arms of a molecule i, with the appropriate orientation 
(fio- ik ,o- u — 1), there is an energy gain J a > 0. 

For J„ — we recover the model of Ref. that 
predicts the singularity-free scenario, and where the HBs 
are uncorrelated, inhibiting the orientational long-range 
order. We study the general case with finite J a , by using 
(i) a mean field approximation and (ii) MC simulations. 



III. THE EQUATION OF STATE APPROACH 



The equation of state of our system is implicitly given 



by 



where 



U ~TS + PV = nj^'' 



U = TL + TLi m 



(5) 



(6) 



is the total internal energy and /i is the chemical poten- 
tial. From Eqs. we rewrite the equation of state 
as 

TS - PV = - ]T 4-(P, a)n inj - £ ^(a)m . (7) 



•:.i> 



Here 



e^.(P )CT ) = e + J'(P)^, CTii 



(8) 



is the effective attraction energy, depending on P and the 
local arm configuration, 



J'(P) = J-PV HB 



(9) 



is the effective HB interaction energy due to the HB 
volume-increase, and 



MiO) = M + Jo S *ki,<rii 
{k,l)i 



(10) 



is the effective local chemical potential depending on the 
local arm configuration. 



A. The mean field approximation 

The mean field approximation consists in assuming a 
linear relation between the number density of liquid cells, 



N 



(11) 



and the density order parameter m G [—1,1], and be- 
tween the number density of arms in the appropriate 
state for a HB, 



(12) 



and the orientational order parameter m a € [0, 1], i.e. 



1 + m 



1 + (q - l)m a 



7 (J 

2 q 
Hence, the molar density p = nN/V is 
1 + m 



(13) 



(14) 



2v + AvhbPhb 

Here = Vo/N, and phb = n 2 p a is the probability of 
forming a HB between two nn molecules, where n 2 is the 
probability of finding two nn molecules and 



p a = n a + (q - 1) 



1 - n a 

<z-i 



\ + {q-l)m 2 a 



(15) 



is the probability of having the facing arms of the two 
molecules in the appropriate orientational state for a HB. 

For T — > oo we expect m a — * 0, hence p a — > 1/q. For 
T — > the finite values of e, J and J CT allow us to assume 
a cooperative effect and an orientational long-range order 
in a preferred state, with m a — > 1, hence p a — > 1. 
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B. The cooperative effect 



C. The Gibbs free energy 



To include the cooperative effect, we consider that each 
arm cr^ interacts with a mean field h generated by the 
other three arms on the same molecule, in addition to the 
effective interaction J'(P) with the facing arm (jji on a nn 
molecule. Since the energy is minimized when the arms 
are in the same orientational state, the system breaks the 
symmetry ordering in the preferred Potts state. Hence, 
we choose h proportional to n a , to J a and to the number 
of arms generating h, i.e. 



h = ZJ a n a 



(16) 



Our results do not depend on this choice for h, and are 
recovered using higher order approximations for h. 
We assume 



{S„ 



(17) 



i. e. that p a [Eq.(tL5|)], for two nn molecules interact- 
ing with the surrounding, coincides with the probability 
(£<Ti -,a-i)h, f° r the facing arms (cr.^ , <jji) of two isolated nn 
molecules, of being in the same orientational state under 
the action of the field h. By definition is 



(So 



= jT E 5 ° ,c exp{[J'(P)S a , T 



Z h 

<7ij ,< 

1 + (g - 1 



8^,i)]/(k B T)} 
2w ma +q-2 



(18) 
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where the right-most-side is the explicit calculation of the 
left-side with partition function Zh- Here the sum is over 
all the configurations of the two variables cry, o~ji, and 
the symbols are 



Next, we write a mean field expression for the molar 
Gibbs free energy 



g = u — Ts + Pv — fi 



(20) 



as a function of the two order parameters m and to*, 
where 

u=^-=-2[en+{Jn + 2,J a )p a ] (21) 

is the molar energy as derived by Eqs. and Eq. (||), 

with the mean field approximations 

riiUj — n 2 , 

(id) 

11,1' , CT = Tl 2 p a , 

(id) 

^ S <r t k,<ra = n Va \ 
i (k,l)i 



(22) 



v — 1/p is the molar volume derived by Eq. 

Sw + S a 



nN 



(23) 



is the molar entropy, 
Sw 



k B N 
Ak B nN 



nlnn + (1 — n) ln(l — n) 



(24) 



. 1 — n* \ f 1 — nZ 
= <ln< + (g-l) [ & I In 
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are the standard mean field expressions for the entropy 
of N variables rc.j and for the entropy of 4nN g-states 
Potts variables for the arms, respectively, and n* is the 
number density n a of HBonded arms [Eq. (pHI)] evaluated 
in m*(T,P). 



U'n 



exp 



exp 



J'(P) 



3J a [l+m a (q-l)] 



E ex P 

7 ij t°"ji 



qk B T 

J'{P)5a 



(19) 



with the Boltzmann constant ks chosen as unitary here- 
after. As expected for p a , also {8 aijlC rn)h — > l/<7 for 
T — > 00 and (5 - y , < r ii )ft — > 1 for T — * 0. 

Numerically we find that the solution of the Eq. ([l7]) is 
m* a {T,P) = for P > P max {T), and m* a {T,P) > for 
P < Pmax{T), where m*(T,P) = corresponds to lack 
of orientational order and m* (T, P) > corresponds to 
the orientational long-range order. In this mean field 
approximation, P max (T) turns out to be well described 
by a decreasing linear function of T. 



D. The mean field results 

By numerically minimizing g(T, P) with respect to 
m and to* with the constraint that to* is solution of 
Eq. jl7|), we find the equilibrium values of m(T,P) and 
m* (T, P). By using Eq. |TJ), we find of p(T, P) at equi- 
librium (Fig. |). 

At high P the mean field theory predicts that p(T) 
increases when T decreases (Fig. 1). At low P, for de- 
creasing T, the theory predicts (i) a discontinuity in p(T), 
corresponding to the liquid-gas first order phase transi- 
tion ending in the liquid-gas critical point C (Fig. |2|), (ii) 
the TMD decreasing for increasing P, (in) a discontinu- 
ity in p(T) at low T, disappearing at lower P. 

The first two predictions are consistent with either the 
singularity-free scenario or the liquid-liquid phase tran- 
sition hypothesis, while the third prediction is consis- 
tent only with the HDL-LDL first order phase transition 
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FIG. 1: The mean field isobaric molar density p as function 
of T for (top to bottom) Pv /e = 1.6, 1.4, 1.3, 1.275, 1.0, 
0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.25, 0.2, 0.15, 0.1, 0, for the 
model with parameters q = 6, J/e = 0.5, J CT /e = 0.05, and 
vhb/vq = 0.5. (a) For 0.1 < Pvo/e < 0.25, by decreasing 
T, p has a discontinuity at high T, from a low value (in the 
gas phase) to a high value (in the liquid phase), then p has a 
maximum followed by a smooth saturation to the finite value 
p HB — 0.5/«o corresponding to the full-HBonded liquid. For 
0.3 < Pvo/e < 0.4, by decreasing T, there is no discontinuity 
in p, but there is a maximum in p and the saturation to pus- 
For 0.5 < Pvo/e < 1.2, by decreasing T, p has a maximum 
and then a discontinuity to phb- For 1.25 < Pvo/e < 1.6, p 
has only a maximum, and, for higher P, p regularly increases, 
by decreasing T. (b) Blowup of the low-T region. Both dis- 
continuities reported show a first order phase transition, each 
ending in a critical point (Fig. 2). 

hypothesized in the latter scenario. In particular, the 
smooth disappearing of the discontinuity at lower P is 
consistent with a phase transition line, with a negative 
slope in the P-T phase diagram, ending in a HDL-LDL 
critical point C (Fig. |). 

IV. THE OFF-LATTICE MODEL 

To show that our mean field predictions are robust, we 
now use a completely different approach based on an off- 
lattice (OL) model representing a system with a homo- 
geneous distribution of molecules in the available volume 
V, that we divide in N equivalent cells of volume V/N. 

As a consequence of the homogeneity of the system, 
for each cell the N degrees of occupancy freedom m are 
set to rii — 1. In analogy with Eq. (g), the total volume 
V is defined as 

V = V OL + N HB v HB . (25) 

Here Nhb and vhb are defined as in Eqs. but, 
differently from the lattice case, the volume , asso- 
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0.5 1 

Temperature T/e 



FIG. 2: The P-T phase diagram for the model with the pa- 
rameters in Fig. 1. The squares and the circles are estimated 
from the discontinuities and the maxima in Fig. 1, respec- 
tively. Since the symmetry between the two phases separated 
by a critical point must be preserved, the high-T discontinuity 
shows a gas-liquid phase transition, while the low-T discon- 
tinuity shows the HDL-LDL phase transition. The lines are 
guides for the eyes. The liquid-gas phase transition line ends 
in the critical point C. The HDL-LDL phase transition line, 
with negative slope, ends in the critical point C". The dashed 
line corresponds to the TMD line. 



dated with the total volume of the cells without HBs, is 
a continuous variable with the constrains V® L > Nvq, 
where vq is the hard-core volume of a molecule. 

Also, following the lattice model, the molecules have 
four arms described by four g-state variables cry, with 
the HB interaction defined by the second term in Eq. (|l|) 
and the intra-molecular interaction by Eq. (|j). These 
interactions are both independent of the distance among 
first-neighbor molecules and depend only on the arms 
orientation aij . 

In this off-lattice model, the average distance between 
two molecules r is a continuous variable, so we replace 
the first term in Eq. (jl]) with 

{oo for r < Rq 

4(^) 12 -(^) 6 ] ^ r>R , W 

where Rq = ^/vo is the hard-core diameter of each 
molecule. In analogy with the Eq. (|l|), we consider this 
off-lattice van der Waals energy independent of the HB 
expansion, so in Uw(r) we use 

r ee y/v° L /N . (27) 
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FIG. 3: MC isobaric density p MC (T) for iV = 10 4 molecules, 
for the off- lattice model with parameters as in Fig. 1. We 
show only the isobars for (bottom to top) Pvo/e = 0.1, 0.2, 
0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95, 0.975, 1, 1.1. (a) The 
qualitative behavior is as described in Fig. 1. (b) Blowup of 
the low-T region. 



A. The Monte Carlo simulation 

We perform MC simulations, in two dimensions p9| , at 
constant N, P, T and variable V (NPT ensemblejwith 
iV G [10 2 ,10 4 ]. The MC dynamics consists in updating 
the variables atj and the variable Vq >l , accepting the new 
state with probability exp[-(AU w +PAV/k B T)} if AU+ 
PAV > 0, or with probability 1 if AU + PAV < 0. Here 
AU = A{Uw + 'Him) and AV are the changes of total 
internal energy and total volume Eq.(p5|), respectively, 
after the update. Our results for the average density 
p MC = N jy^ averagec i n 6x 10 5 MC steps after 1.2 x 10 5 
MC steps of thermalization at each T, are qualitatively 
consistent with the mean field prediction (Fig. ||). 

By MC simulations, we find (i) for J a = 0, no liquid- 
liquid phase transition and the TMD line, consistent with 
mean field in Rcf. p2] j; (ii) for J a ^ J, for any P in the 
liquid phase, at T below the TMD line, a discontinuity 
in density, suggesting a first order phase transition along 
a line with a negative slope in the P-T phase diagram, 
consistent with mean field for J a — > oo |2?j, ^8); (iii) for 
< J a < J, a phase transition line ending in a criti- 
cal point C (Fig. [|), and occurring at increasing P and 
decreasing T for decreasing J a . 

To verify that the jumps found in the MC density are 
marking a first order phase transition, instead of a nar- 
row continuous phase transition, we study the isothermal 
compressibility, 




I , 1 , , , 1 

1 2 

Temperature T/e 

FIG. 4: P-T phase diagram calculated by MC simulations, 
for the off-lattice model as in Fig. 3. Squares are the N — > oo 
estimates for the points on the coexisting lines. Full circles are 
the critical points C and C". Points on the TMD line (open 
circles) are estimated from the N — > oo extrapolation of the 
maxima of p MC (T, N). Dashed lines indicate the position 
of K^ ax emanating from C and C' . Full lines are guides 
for the eyes. Where not shown, errors are smaller then the 
symbol size. Inset: blowup of the HDL-LDL phase transition 
region. The full line is given by the empiric expression P ma x = 
Pa - [a/ (T - T)] where P v /e = 1.087, av /e 2 = 0.006 and 
T /e = 0.07. 



Its maximum K™ ax , where 

(dK T /dV)\ T = 
(d 2 K T /dV 2 )\ T > , (29) 

increases linearly with the number of particles N at a 
first order phase transition |?0), while at a second order 
phase transition K™ ax is proportional to both N and the 
fluctuation of the ordering parameter, scaling as a power 
of N pl| . Therefore, the finite size scaling analysis on 
K™ ax (N) allows us to discriminate between a continuous 
and a discontinuous phase transition. 

To estimate the maximum with a great precision, we 
use a continuous T algorithm, the histogram reweighting 
method p2[ . By checking the minimum T and the maxi- 
mum P where the behavior of K™ ax (N) fails to be linear, 
we estimate the critical point at T"c'/ e = 0.045 ± 0.005 
and Pc'Vo/e = 0.841 ± 0.042 for N -> oo (inset Fig.4). 

Next, we obtain the coexisting lines by extrapolating to 
N -> oo the values P(T,N) corresponding to K^ ax (N) 
at fixed T, both for T < T c and T < T c > . Furthermore, 
our results are consistent with the necessary condition 
that the K™ ax (T) line emanates from the critical point 
H |§, both for T > T C ' and for T > T c (Fig.4). 
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V. DISCUSSION AND CONCLUSIONS 

The mean field and the MC results for our water-like 
fluid model for a finite intra-molecular interaction J a 
show qualitatively the same phase diagram. In both ap- 
proaches the TMD line decreases with increasing P, con- 
sistent with the experiments 0]. In particular, both ap- 
proaches predict a first order phase transition in the liq- 
uid phase, occurring at low T and at a pressure P m ax{T) 
decreasing for increasing T. For J CT < J we find that the 
first order phase transition line ends in a critical point, 
separating, by necessity, two phases with the same sym- 
metry, in this case two liquids, HDL and LDL, being the 
critical point in the liquid phase. Three considerations 
are in order here. 

The first consideration is related to the comparison 
with the result of Ref. p2|, h ere recovered for J a = 0. 
The Sastry et al. model p2|], upon HB formation, ac- 
counts for (i) inter-molecular orientational correlation 
[Eq.(l)]; (ii) local expansion with lowering temperature 
[Eq.(2)]; (Hi) anticorrelation between V and S, because 
the formation of HBs decreases the number of possible 
orientational configurations for the system, hence the en- 
tropy S decreases for increasing Nhb, i-e. for increasing 
V . This is expected in a system with a density anomaly, 
because [dV/dT) P < implies [dS/dV) T < 0. Fi- 
nally, the Sastry et al. model assumes the arms of a 
molecule completely independent (J a — 0), and predicts 
the singularity- free scenario. 

We tested, by preliminary MC calculations, that for 
J a — > the HDL-LDL critical point C moves to lower T 
and to P m ax{T = 0). In particular, the phase transition 
disappears at T = for J a = 0, while the effect of the 
decreasing J a on the location of the TMD line is weak. 
Although these results require a longer analysis, beyond 
the scope of the present work, they show that the pre- 
dictions of Ref. [^] are recovered in the limit J a = 0, 
confirming the validity of our MC approach |j4j. We, 
therefore, conclude that in this model the presence of a 
finite J iy is responsible for the appearance of the first or- 
der phase transition, with a possible HDL-LDL critical 
point C. 

The second consideration is related to the possibility 
that this low-T HDL-LDL phase transition is pre-empted 
by inevitable freezing in real water. Recent analysis of 
the realistic model for water ST2 [H| suggests that the 
HDL-LDL critical point may occur above the glass tem- 
perature T g , though as yet still outside the easily acces- 
sible experimental range. 

We compare our MC results with data for real water. 
From the location of the liquid-gas critical point (Tc, Pc) 
and the TMD line at ambient pressure (T*,P*) (Table |), 
we find the ratios P* /Pc and T*/T c in real water. By 
assuming that the same P* /Pc holds in our MC case, we 
calculate the corresponding P* ~ 2.76 10~ 4 e/vo in our 
model and then we estimate the T* corresponding to P* 
from the TMD line in the MC phase diagram (Fig. ||). 
In this way we find a ratio T* /Tc from the MC results 



that is consistent with the real water data, suggesting the 
validity of our assumption on P*/Pc- 

Therefore, we use the same kind of assumption also 
to estimate the glass temperature T g for our phase di- 
agram. In particular, from real water data we obtain 
the ratio T g /Tc at ambient pressure and, assuming that 
it holds also for our model, we estimate T g /e ~ 0.42 at 
Pvg/e ~ 2.76 10~ 4 for the MC phase diagram. Hence, 
for our model with the parameters chosen in this pa- 
per, is T g > Tc, i.e. the HDL-LDL critical temperature 
at Pc — 0.841e/i>o is below the glass temperature at 
P* ~ 2.76 10~ 4 e/t! . From the study of the phase dia- 
gram of real water Q is reasonable to expect that T g (P) 
decreases for increasing P, therefore our analysis does 
not exclude that Tc is above T g (Pc). However, by con- 
sidering a very large J a , such that the HDL-LDL critical 
pressure is ~ 2.76 10~ 4 e/t>o, we can compare Tc and T g 
at the same pressure P*. Our preliminary results show 
that the HDL-LDL critical temperature is in this case 
very close to T g . 

As a consequence of this analysis, our model supports 
the possibility that the HDL-LDL critical point is lo- 
cated deep into the supercooled region, below or close 
to the glass temperature, depending on the value of J a . 
Therefore, the liquid-liquid phase transition could be pre- 
empted, if our model is representative of the thermody- 
namic properties of the real system. This results is anal- 
ogous to what has been proposed for silica p5| , another 
liquid with density anomaly, suggesting that the present 
model could provide a general theoretical framework for 
anomalous molecular liquids. 

The third consideration is about the role of the tetra- 
hedrality in determining the properties of anomalous liq- 
uids. For the present model we do not consider a tetra- 
hedral geometry in the two-dimensional MC approach, 
and the geometry is not explicitly defined in the mean 
field approach. Nevertheless, our results are consistent 
with the experimentally accessible phase diagram of real 
water, suggesting that the tetrahedral network is not an 
essential feature for the anomalous behavior of water- like 
liquids. 

This conclusion is consistent with what has been ob- 
served by Angell in Ref. and is well described by a 
general cooperative model Q with a generic drive to 
phase separate the excitations into distinct regions of 
space (clustering). In our model the drive is given by 
the intra-molecular interaction, that mimics the geomet- 
rical drive in tetrahedral liquids even if is not necessarily 
limited to the tetrahedral case. 

In conclusion, we studied the effect of an intra- 
molecular interaction J a in a model for anomalous molec- 
ular liquids with a mean field approach, valid for J a > 0, 
and with an off-lattice MC simulation. For J„ > we 
found a HDL-LDL phase transition while our MC re- 
sults confirm that for J a = the singularity-free sce- 
nario holds p^ ). Hence, the two interpretations originate 
from the same mechanism with a different hypothesis on 
the intra-molecular interaction; the latter is strictly valid 
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TABLE I: Characteristic temperatures and pressures for real water and for the present model, for the gas-liquid critical point 
(Tc,Pc), the TMD at ambient pressure (T*,P*) and the glass temperature T g at ambient pressure. The ratios P*/Pc and 
Tg/Tc are not available for the model. We assume that the corresponding H2O values are valid also for the model, and we 
use these values to estimate P* and T g , respectively, for the model. Temperatures are measured in K for H2O and in e for the 
model. Pressures are measured in MPa for H2O and in e/vo for the model. 





Pc 


T C 




pyp c 


T* 


T*/T C 


T a T g /Tc 


H 2 


22.064 


647.14 


0.10133 Q 


4.6 10~ 4 


277 |5) 


0.428 


136 (36) 0.21 


model 


0.6 ±0.1 


2.0 ±0.2 


~ 2.76 10 -4 




0.7 ±0.1 


0.35 ± 0.06 


~ 0.42 



only for J a — 0. Within this framework, the most rea- 
sonable scenario for water includes a HDL-LDL phase 
transition, probably hindered by inevitable freezing. Our 
results suggest also that the tetrahedrality is not essen- 
tial to understand the anomalous behavior in water-like 
liquids. 
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